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Abstract. High angular resolution H 2 CO 218 GHz line observations have been carried out toward the low-mass protostars 
IRAS 16293-2422 andL1448-C using the Owens Valley Millimeter Array at ~2" resolution. Simultaneous 1.37 mm continuum 
data reveal extended emission which is compared with that predicted by model envelopes constrained from single-dish data. 
For L1448-C the model density structure works well down to the 400 AU scale to which the interferometer is sensitive. For 
IRAS 16293-2422, a known proto-binary object, the interferometer observations indicate that the binary has cleared much 
. of the material in the inner part of the envelope, out to the binary separation of ~800 AU. For both sources there is excess 

' unresolved compact emission centered on the sources, most likely due to accretion disks <200AU in size with masses of 

>0.02 M Q (L1448-C) and >0.1 M (IRAS 16293-2422). The H 2 CO data for both sources are dominated by emission from 
' gas close to the positions of the continuum peaks. The morphology and velocity structure of the H2CO array data have been 

used to investigate whether the abundance enhancements inferred from single-dish modelling are due to thermal evaporation of 
ices or due to liberation of the ice mantles by shocks in the inner envelope. For IRAS 16293-2422 the H 2 CO interferometer 
observations indicate the presence of large scale rotation roughly perpendicular to the large scale CO outflow. The H 2 CO 
distribution differs from that of C ls O, with C 18 emission peaking near MM1 and H 2 CO stronger near MM2. For L1448-C, 
the region of enhanced H 2 CO emission extends over a much larger scale > 1" than the radius of 50 - 100 K (076 - 07 15) where 
thermal evaporation can occur. The red-blue asymmetry of the emission is consistent with the outflow; however the velocities 
JL , are significantly lower. The H 2 CO 3 22 - 2 21 /3o3 - 2 02 flux ratio derived from the interferometer data is significantly higher 

than that found from single-dish observations for both objects, suggesting that the compact emission arises from warmer gas. 
Detailed radiative transfer modeling shows, however, that the ratio is affected by abundance gradients and optical depth in the 
3o3 _ 2 02 line. It is concluded that a constant H 2 CO abundance throughout the envelope cannot fit the interferometer data of 
the two H 2 CO lines simultaneously on the longest and shortest baselines. A scenario in which the H 2 CO abundance drops in 
the cold dense part of the envelope where CO is frozen out but is undepleted in the outermost region provides good fits to the 
single-dish and interferometer data on short baselines for both sources. Emission on the longer baselines is best reproduced 
5^ ' if the H 2 CO abundance is increased by about an order of magnitude from ~10~ 10 to ~10~ 9 in the inner parts of the envelope 

due to thermal evaporation when the temperature exceeds ~50 K. The presence of additional H 2 CO abundance jumps in the 
innermost hot core region or in the disk cannot be firmly established, however, with the present sensitivity and resolution. Other 
scenarios, including weak outflow-envelope interactions and photon heating of the envelope, are discussed and predictions for 
future generation interferometers are presented, illustrating their potential in distinguishing these competing scenarios. 

Key words, astrochemistry - stars: formation - circumstellar matter - stars: individual: IRAS 16293-2422, L1448-C - ISM: 
abundances 



1. Introduction 

Recent observational studies have shown that the inner (< 
few hundred AU) envelopes of low-mass protostars are dense 
(>10 6 cm 4 ) and warm (>80 K) (Blake et al. 1994; Ceccarelli 
et al. 2000a; J0rgensen et al. 2002; Schoier et al. 2002; Shirley 
et al. 2002), as expected from scaling of high-mass protostars 
(Ceccarelli et al. 1996; Ivezic & Elitzur 1997). In high-mass 
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objects, these warm and dense regions are known to have a 
rich chemistry with high abundances of organic molecules due 
to the thermal evaporation of ices (e.g., Blake et al. 1987; 
Charnley et al. 1992). Detailed modeling of multi-transition 
single-dish lines toward the deeply embedded low-mass pro- 
tostar IRAS 16293-2422 has demonstrated that similar en- 
hancements of molecules like H2CO and CH3OH can occur 
for low-mass objects (van Dishoeck et al. 1995; Ceccarelli 
et al. 2000b; Schoier et al. 2002). Recently, Maret et al. (2004) 
have suggested that this is a common phenomenon in low-mass 
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protostars. The location at which this enhancement occurs is 
consistent with the radius at which ices are expected to ther- 
mally evaporate off the grains (T > 90 K). Moreover, large 
organic molecules have recently been detected toward IRAS 
16293-2422 (Cazaux et al. 2003), showing that low-mass hot 
cores may have a similar chemical complexity as the high-mass 
counterparts in spite of their much shorter timescales (Schoier 
et al. 2002). 

Alternatively, shocks due to the interaction of the outflow 
with the inner envelope can liberate grain mantle material over 
a larger area than can thermal heating. Additionally, the bipolar 
outflow will excavate a biconical cavity in the envelope through 
which UV- and X-ray photons can escape. The back scattering 
of such photons into the envelope by low-density dust in the 
cavity can significantly heat the envelope surrounding the cav- 
ity (e.g., Spaans et al. 1995). This would produce regions of 
warm gas (~100 K) in the envelope on much larger scales than 
otherwise possible. High angular resolution observations are 
needed to pinpoint the origin of the abundance enhancements 
and distinguish between these various scenarios. 

We present here observations of H2CO toward two low- 
mass protostars, IRAS 16293-2422 and L1448-C (also known 
as L1448-mm), at 218 GHz (1.4 mm) using the Owens Valley 
Radio Observatory (OVRO) Millimeter Array at ~2" reso- 
lution. The frequency setting includes the 3o3 — > 2o2 and 
322 — » 221 H2CO lines, whose ratio is a measure of the gas 
temperature of the circumstellar material. Both IRAS 16293- 
2422 and L1448-C are deeply embedded class protostars 
(Andre et al. 1993) which drive large scale (~arcmin) bipo- 
lar outflows (Walker et al. 1988; Mizuno et al. 1990; Bachiller 
et al. 1990, 1991; Stark et al. 2004). For other low-mass ob- 
jects, molecules such as SiO are clearly associated with the 
outflow (e.g., L1448:Guilloteauetal. 1992; NGC 1333IRAS4: 
Blake 1995), whereas optically thick lines from other species 
such as HCO + and HCN are found to 'coat' the outflow walls 
(e.g., B5 IRS1: Langeret al. 1996; L1527 and Serpens SMM1: 
Hogerheijde et al. 1997, 1999). The extent of this emission 
can be larger than 10", which should be readily distinguishable 
from the ~1" hot inner envelope with current interferometers. 

Previous millimeter aperture synthesis observations of 
IRAS 16293-2422 have revealed two compact components co- 
incident with radio continuum emission, indicative of a proto- 
binary source (Mundy et al. 1990, 1992). The line emission of 
10 molecular species at ~5" resolution reveals that there is a 
red-blue asymmetry indicative of rotation perpendicular to the 
outflow direction (Schoier et al., in prep.). The morphology of 
the emission picked up by the interferometer suggests that it 
may be produced in regions of compressed gas as a result of in- 
teraction between the outflow and the envelope. Previous data 
on L1448-C show a compact continuum source at millimetre 
wavelengths and that SiO is a good tracer of the large velocity 
outflow associated with this source (Guilloteau et al. 1992). 

Since most of the extended emission is resolved out by 
the interferometer, a good physical and chemical model of 
the envelope is a prerequisite for a thorough interpretation of 
the aperture synthesis data. In recent years, much progress 
has been made in obtaining reliable descriptions of the den- 
sity and temperature structures in the dusty envelopes around 



Table 1. Source and envelope parameters for IRAS 16293-2422 and 
L1448-C 





TP AC 1 fi.lQV 1 


L1448* 


Distance, D (pc) 


160 


220 


Luminosity, L (Lq) 


27 


5 


Inner radius, r, (AU) 


32.1 


9.0 


Outer radius, r e (10 4 AU) 


0.80 


0.81 


Density power law index, a 


1.7 


1.4 


Density at 1000 AU, « (H 2 ) (10 6 cnr 3 ) 


6.7 


0.75 


Column density c , N(H 2 ) (10 24 cnr 2 ) 


1.6 


0.17 


Envelope mass r , M env (Mq) 


5.4 


0.93 



" From Schoier et al. (2002) 
* From l0rgensen et al. (2002) 
c Within the outer radius r e 



young stellar objects, based on thermal continuum emission 
(Chandler & Richer 2000; Hogerheijde & Sandell 2000; Motte 
& Andre 2001; J0rgensen et al. 2002; Schoier et al. 2002; 
Shirley et al. 2002). The physical structures of IRAS 16293- 
2422 and L1448-C have recently been derived from single- 
dish continuum observations, with the results summarized in 
Table 1 (J0rgensen et al. 2002; Schoier et al. 2002). 

In §3, we first test the validity of these envelope models 
at the small scales sampled by the continuum interferometer 
data. In §4, the H2CO results are presented and analyzed. For 
IRAS 16293-2422, C ls O observations are also available. This 
is followed by a discussion on the origin of H2CO and estimates 
for what future generation telescopes might reveal in §5 and by 
conclusions in §6. 

2. Observations and data reduction 

2. 1 . Interferometer data 

The two protostars IRAS 16293-2422 (a 2 ooo = 16 h 32 m 22 s .8, 
£2000 = -24°28'33'/0) and L1448-C (a 2 ooo = 3 h 25 m 38 s .8, 
£2000 = 30°44'05'.'0) were observed with the Owens 
Valley Radio Observatory (OVRO) Millimeter Array 1 between 
September 2000 and March 2002. The H 2 CO 3 03 -» 2 02 and 
322 -> 2 2 i line emission at 218.222 and 218.475 GHz, respec- 
tively, was obtained simultaneously with the continuum emis- 
sion at 1.37 mm. IRAS 16293-2422 was observed in the L 
and E configurations, while L1448-C was observed in the C, 
L and H configurations, corresponding to projected baselines 
of 8 - 80 and 8-120 kA, respectively. The complex gains 
were calibrated by regular observations of the quasars NRAO 
530 and 1622-253 for IRAS 16293-2422 and 0234+285 for 
L1448-C, while flux calibration was done using observations 
of Uranus and Neptune for each track, both using the MMA 
package developed for OVRO data by Scoville et al. (1993). 

1 Research with the Owens Valley Millimeter Array, operated by 
California Institute of Technology, is supported by NSF grant AST 
99-81546. 
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The subsequent data-reduction and analysis was performed us- 
ing MIRIAD (Sault et al. 1995). 

Further reduction of the data was carried out using the stan- 
dard approach by flagging clearly deviating phases and ampli- 
tudes. The continuum data were self-calibrated and the result- 
ing phase corrections were applied to the spectral line data, 
optimizing the signal-to-noise. The natural-weighted contin- 
uum observations for IRAS 16293-2422 and L1448-C have 
typical lcr noise levels better than 20 and 3mJybeam _1 with 
beam sizes of 3'.'9xr.'9 and 2'.'6x2'.'3, respectively. The rela- 
tively high noise levels for IRAS 16293-2422 reflect the low 
elevation at which this source is observable from Owens Valley, 
which both increases the system temperatures and decreases 
the time available per transit. The data were then CLEANed 
(Hogbom 1974) down to the 2<x noise level. 

For IRAS 16293-2422 additional archival C ls O 7 = 2^1 
line data obtained in 1993 using the OVRO array are presented, 
which were reduced in the same way as described above. 

2.2. Single-dish data 

It is well-known that the millimeter aperture synthesis obser- 
vations lack sensitivity to extended emission due to discrete 
sampling in the (u, v) plane and, in particular, missing short - 
spacings. In order to quantify this single-dish observations 
were performed using the James Clerk Maxwell Telescope 
(JCMT) 2 . The continuum data are taken largely from the JCMT 
archive 3 and have been presented in Schoier et al. (2002) and 
J0rgensen et al. (2002). For H2CO, a 25-point grid centered on 
the adopted source position and sampled at 10" spacing was 
obtained in September 2002 for IRAS 16293-2422, with both 
H2CO 218 GHz lines covered in a single spectral setting. The 
observations were obtained in a beam-switching mode using a 
180" chop throw. The data were calibrated using the chopper- 
wheel method and the resulting antenna temperature was con- 
verted into main-beam brightness temperature, r m b, using the 
main-beam efficiency 77 m b = 0.69. For L1448-C, a single spec- 
trum at the source position was taken that includes both transi- 
tions. 

3. Continuum emission: disk and envelope 
structure 

3.1. L1448-C 

In Fig. 1, the 221.7GHz (1.37mm) continuum emission to- 
ward L1448-C is presented. Only a single compact compo- 
nent is seen, with faint extended emission. The total continuum 
flux density at 1.37 mm observed with OVRO is 0.32 Jy, only 
35% of the flux observed by Motte & Andre (2001) (0.9 Jy) at 
1 .3 mm using the IRAM 30 m telescope. The compact compo- 
nent, located at (—076, 072) from the pointing centre, has been 

2 The JCMT is operated by the Joint Astronomy Centre in Hilo, 
Hawaii on behalf of the Particle Physics and Astronomy Research 
Council in the United Kingdom, the National Research Council of 
Canada and the Netherlands Organization for Scientific Research. 

3 http://www. jach.hawaii . edu/JACpublic/JCMT/ 
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Fig. 1. OVRO interferometer maps of the 1.37 mm continuum emis- 
sion toward L1448-C. Contours start at 9 mJybeanT 1 and the first 
contour corresponds to the '2cr'-level as estimated from the maps. 
Each successive contour is a multiple of two of the preceding value. 
The peak value is 0. 13 Jy beam -1 and the synthesized beam in this and 
subsequent images is depicted by the filled ellipse in the lower-right 
corner. 





0.4 










: L1448-C (1.37 


mm) : 










'" a 




Envelope with a=1.4 




0) 


0.3 


Envelope with 01=1.4 + point 




p 




* Envelope with ot=l.& + point 


source 










amp] 


0.2 






Llity 




: ^*>>- r -. 1 

: • ^^i^ui-- 






0.1 






\ISI 


; + 






0.0 







20 40 60 80 

uv distance [kX] 

Fig. 2. Visibility amplitudes of the observed 1.37 mm continuum emis- 
sion obtained at OVRO towards L1448-C as functions of the projected 
baseline length, binned to 5 kA, from the phase center, taken to be at 
(-076,072). The observations are plotted as filled symbols with lcr 
error bars. The dotted histogram represents the zero-expectation level. 
Also shown are predictions based on a realistic physical model for 
the source (J0rgensen et al. 2002), with the same (u, v) sampling as 
the observations (see text for details). Unresolved compact emission, 
presumably from a circumstellar disk, must be added to that from the 
envelope to produce an acceptable fit. 

fitted with a Gaussian in the (u, v) plane, resulting in an esti- 
mated size of 170x076 and an upper limit to the diameter of 
~170 AU for the adopted distance of 220 pc. 

J0rgensen et al. (2002) determined the actual tempera- 
ture and density distribution of the circumstellar envelope of 
L1448-C from detailed modeling of the observed continuum 
emission (see Table 1). In addition to the spectral energy dis- 
tribution (SED), resolved images at 450 and 850 i/m obtained 
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with the SCUBA bolometer array at the JCMT were used to 
constrain the large scale envelope structure. The interferome- 
ter data constrain the envelope structure at smaller scales (~2") 
than the JCMT single-dish data (~10-20")- In order to inves- 
tigate whether this envelope model can be reconciled with the 
flux picked up by the interferometer, the same (u, v) sampling 
was applied to the predicted brightness distribution at 1.37 mm 
from the model envelope. 

Fig. 2 compares the observed visibility amplitudes of 
L1448-C with the model predictions. First, it is clear that the 
envelope alone cannot reproduce the data, but that an unre- 
solved compact source, presumably the disk, needs to be added. 
The combination of this point source and the best fit envelope 
model of J0rgensen et al. (2002) with a density structure falling 
off as r~ a , where a — 1.4, dramatically improves the fit to the 
visibilities in the (u, v) plane. A slightly steeper density struc- 
ture of a — 1 .6 is preferred by the interferometer data. Given 
the mutual uncertainties of +0.2 in a these results are still in 
good agreement with each other and with those of Shirley et al. 
(2002), who found a slope a — 1 .7 in their analysis for a some- 
what larger (r e = 45000 AU) envelope. The remaining point 
source flux for L1448-C is estimated to be 100 and 75 mJy 
when a is taken to be 1.4 and 1.6, respectively. This is ~ 10% 
of the total single-dish source flux. The impact of changing dif- 
ferent envelope parameters such as density slope and location 
of inner envelope radius was tested in J0rgensen et al. (2004a), 
for the embedded low-mass protostar NGC 1 333-IRAS2A, and 
found to lead to an uncertainty of about +25% on the derived 
point source flux. 

The point source flux estimated here at 1.37 mm agrees 
very well with the spectral index of 1.84 + 0.08 found from 
other mm and cm wavelength observations (Curiel et al. 1990; 
Guilloteau et al. 1992; Looney et al. 2000; Reipurth et al. 2002) 
as shown in Fig. 3. The quoted standard deviation is based 
on assuming a 20% uncertainty in all the fluxes. Similarly, 
J0rgensen et al. (2004a) found a spectral index of 1.9 for 
the point source associated with the low-mass protostar NGC 
1333-IRAS2A. The spectral index is consistent with optically 
thick thermal emission and its favoured origin is that from an 
unresolved accretion disk. 

Assuming the point source emission to be thermal the mass 
of the compact region can be estimated from 



M 



KyBy(Td) 

where F v is the flux, *F is the gas-to-dust ratio (assumed to be 
equal to 100), D is the distance, k v is the dust opacity, B v is 
the Planck function at a characteristic dust temperature and 
t v is the optical depth. The adopted dust opacity at 1 .37 mm, 
0.8 cm 2 g _1 , is extrapolated from the opacities presented by 
Ossenkopf & Henning (1994) for grains with thin ice man- 
tles. These opacities were used also in the radiative transfer 
analysis of the envelope. For a dust temperature in the range 
100-40 K, the estimated disk mass in the optically thin limit 
is 0.016 - 0.042 M Q when a point source flux of 75 mJy is 
used. This should be treated as a lower limit since the emission 
is likely to be optically thick at 1.37 mm, as suggested by the 
spectral index. 
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Fig. 3. Continuum flux observations of the compact emission toward 
L1448-C (squares with error bars). The solid line shows a fit to the 
data using F oc i/ 5 , where /3 = 1.84 ± 0.08 is the spectral index. 
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Fig. 4. OVRO interferometer maps of the continuum emission at 
1.37 mm towards IRAS 16293-2422. Contours start at 60 mJy beairT 1 
and the first contour corresponds to the '2cr'-level as estimated from 
the maps. Each successive contour is a multiple of two of the preced- 
ing value. The emission peaks at 0.82 and 1.1 Jy beairr' at MM1 and 
MM2, respectively. 



It is difficult to estimate accurate disk masses for deeply 
embedded sources since it involves a good knowledge about 
the envelope structure, in addition to, e.g., the disk tempera- 
ture. For more evolved protostars where the confusion with the 
envelope is less problematic, disk masses of ~0.01 - 0.08 M 
are derived (e.g., Looney et al. 2000; Mundy et al. 2000), com- 
parable to the values found here. 
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Fig. 5. Visibility amplitudes of the observed 1.37 mm emission ob- 
tained at OVRO towards IRAS 16293-2422 as functions of the pro- 
jected baseline length, binned to 2 Li, from the phase center, taken 
to be at (0,0). The observations are plotted as filled symbols with l<x 
error bars. The dotted histogram represents the zero-expectation level. 
Also shown are predictions based on a realistic physical model for 
the source (Schoier et al. 2002), with the same (u, v) sampling as the 
observations (see text for details). Unresolved compact emission, pre- 
sumably from two circumstellar disks, needs to be added to that ema- 
nating from the envelope in order to obtain an acceptable fit. A model 
envelope with a cavity (solid line), in addition to the unresolved com- 
pact emission, is shown to best reproduce the observations. 



3.2. IRAS 16293-2422 

For the proto-binary object IRAS 16293-2422, two unresolved 
continuum sources are detected separated by approximately 5" 
(Fig. 4). The total observed continuum flux density at 1.37 mm 
is about 3.5 Jy. This is roughly 50% of the flux obtained 
from mapping with single-dish telescopes (Walker et al. 1990; 
Andre & Montmerle 1994), indicating that the interferometer 
resolves out some of the emission. The positions of the contin- 
uum sources [(2'.'0, —2'.'9) and (— 1'/6, Of.' 5)] are consistent with 
the two 3 mm sources MM1 (southeast) and MM2 (northwest) 
found by Mundy et al. (1992). At the distance of IRAS 16293- 
2422 (160 pc) the projected separation of the sources is about 
800AU. 

Schoier et al. (2002) have modeled the circumbinary en- 
velope of IRAS 16293-2422 in detail based on SCUBA im- 
ages and the measured SED. Fig. 5 shows that this envelope 
alone cannot fit the compact sources. The addition of two com- 
pact sources, at the locations of MM1 and MM2, to the best- 
fit envelope model of Schoier et al. (2002) produces the cor- 
rect amount of flux at the longest baselines and the smaller 
baselines, but now provides too much emission at intermedi- 
ate baselines (10 - 30 kA, ~10"), i.e., at scales of the binary 
separation. 

It is found that an envelope model which is void of ma- 
terial on scales smaller than the binary separation best repro- 
duces the observed visibilities, in combination with the com- 
pact emission. For this 'cavity' model a slightly steeper density 
profile is obtained, a - 1.9, from re-analyzing the SCUBA 
images and the SED. Also, the temperature is higher within 
r s 1.5 x 10 16 cm (1000 AU) compared to the standard enve- 
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Fig. 6. Continuum flux observations of the compact emission towards 
IRAS 16293-2422 (squares with error bars). The solid lines show fits 
to the data using F oc where [} is the spectral index. For MM1 a 
combination of two powerlaws were used. 



lope, although the temperature never exceeds 80 K in the cavity 
model. 

Theory has shown that an embedded binary system will un- 
dergo tidal truncation and gradually clear its immediate envi- 
ronment due to transfer of angular momentum from the binary 
to the disk. Thus, an inner gap or cavity with very low den- 
sity is produced (e.g., Bate & Bonnell 1997; Giinther & Kley 
2002). Two binary sources in the T Tauri stage have been im- 
aged in great detail; GG Tau and UY Aur (e.g., Dutrey et al. 
1994; Duvert et al. 1998). Wood et al. (1999) estimate that GG 
Tau has cleared its inner 200 AU radius of material and that the 
bulk of material is located in a circumbinary ring of thickness 
600 AU. IRAS 16293-2422 could possibly be a 'GG Tau in the 
making'. 
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Fig. 7. OVRO interferometer maps of H2CO emission (contours) over- 
layed on the 1.37 mm continuum emission (greyscale) for L1448- 
C. The H 2 CO emission has been separated into a red (dashed lines) 
and a blue (solid lines) part (see text for details). Contours start at 
at 0.2 Jy beam -1 kms~' (2cr) and each succesive contour is a multiple 
of this value. Also indicated are the directions of the large scale CO 
outflow. 



The emission from the two unresolved components is esti- 
mated to contribute ~25% (1.8 Jy) to the total flux at 1.37 mm. 
In Fig. 6 flux estimates for the compact components around 
MM1 and MM2 are compared to those at cm to mm wave- 
lengths (Wootten 1989; Estalella et al. 1991; Mundy et al. 
1992; Looney et al. 2000). The emission from MM2 is well fit- 
ted over the entire region using a spectral index of 2.35 + 0.06, 
consistent with thermal emission from an optically thick disk. 
For MM1 a combination of two power laws provides the best 
fit. At shorter wavelengths a spectral index of 2.44 + 0.16 is 
observed, presumably thermal disk emission. At longer wave- 
lengths a much lower index of 0.15 + 0.15 is found consis- 
tent with free-free emission from an ionized stellar wind or jet. 
MM1 is indeed thought to be the source driving the large scale 
outflow associated with IRAS 16293-2422. In comparison no 
active outflow is known for MM2. However, it has been sug- 
gested that MM2 is responsible for a fossilised flow in the E- 
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Fig. 8. Comparison between the HiCO line emission towards L1448- 
C, at the source position, from JCMT single-dish observations (line di- 
agram) and OVRO interferometric observations (histogram) restored 
with the JCMT beam (22"). The OVRO spectrum has been scaled in 
order to account for the flux seen in the JCMT spectrum. 



W direction »10" north of MM1 (see Stark et al. 2004, and 
references therein). 

Gaussian fits to the sizes of these disks in the (u, v) plane 
provide upper limits of »250 AU in diameter for MM1 and 
MM2. Using Eq. 1 in the optically thin limit gives estimates 
of the disk masses for MM1 and MM2 of 0.09 - 0.24 M and 
0.12 - 0.33 M , respectively, again assuming the characteristic 
dust temperature to be 100 - 40 K. Since the spectral indices 
indicate that the emission is optically thick in both cases, these 
masses should be treated as lower limits. 
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4. H 2 CO emission: morphology and abundance 
structure 

4.1. L1448-C 

The maps of the H2CO 3o3 — » 2o2 and 322 — > 221 emis- 
sion toward L1448-C are shown in Fig. 7, separated into blue 
(4-5 km s -1 ) and red (5-6 km s~') components. The 3o3 -» 2 2 
emission appears to be slightly resolved with an extension 
to the south along the direction of the outflow. The H2CO 
322 — > 221 is detected only at the source position. The velocity 
structure hints that the emission is related to the known large 
scale outflow, although the velocities are significantly lower 
than the high velocity (typically 30 - 60kms~ 1 ) outflow seen 
in CO and SiO. 

As for the continuum data, care has to be taken when inter- 
preting the line interferometer maps due to the low sensitivity 
to weak large scale emission. A direct comparison between the 
single-dish spectrum and that obtained from the interferome- 
ter observations restored with the single-dish beam is shown 
in Fig. 8. The interferometer picks up only ~ 10 - 20% of 
the single-dish flux, suggesting that the extended cold mate- 
rial is resolved out by the interferometer and that a hotter, more 
compact, component is predominantly picked up. Within the 
considerable noise and coarsened spectral resolution, the line 
profiles are consistent with those obtained at the JCMT 

The 322 — » 2 2 i/3 ( )3 — » 2 2 line ratio is sensitive to tempera- 
ture (e.g., van Dishoecket al. 1993, Mangum & Wootten 1993), 
especially in the regime of 50 - 200 K. For L1448-C, the inter- 
ferometer data give a ratio of 0.68 ± 0.39 indicating the pres- 
ence of hot gas with T > 70 K. For comparison, the single-dish 
line ratio is 0.12 + 0.04 (Maret et al. 2004), corresponding to 
T w 20 - 30 K. Optical depth effects and abundance variations 
with radius (i.e., temperature) can affect this ratio, however, so 
that more detailed radiative transfer modeling is needed for a 
proper interpretation. 

Just as for the continuum data, the analysis of the H2CO 
interferometer data requires a detailed model of emission from 
the more extended envelope as a starting point. Such a model 
has been presented by Maret et al. (2004) based on multi- 
line single-dish observations. Those data have been re-analyzed 
in this work using the Monte Carlo radiative transfer method 
and molecular data adopted in Schoier et al. (2002). The den- 
sity and temperature structures are taken from J0rgensen et al. 
(2002) (see also §3.1) assuming the gas temperature to be cou- 
pled to that of the dust. The lines are assumed to be broad- 
ened by turbulent motions in addition to thermal line broaden- 
ing. The adopted value of the turbulent velocity is 0.7 kms -1 
(J0rgensen et al. 2002). Considering only the para-FUCO data, 
a good fit (x 2 ed = 0.8) is obtained using a constant para-F^CO 
abundance of 6 x 10" 10 throughout the envelope. A similarly 
good fit (x^ ed = 1.0) can be made to the ortho-H^CO single-dish 
data using an abundance of 9 x 10~ 10 . The uncertainty of these 
abundance estimates is approximately 20% within the adopted 
model. The abundance of H2CO is ~50% larger than that de- 
rived for the outer envelope of IRAS 16293-2422 (Schoier 
et al. 2002). 



As shown by Maret et al. (2004), a different interpretation 
is possible within the same physical model if the ortho-to-para 
ratio is forced to be equal to 3 and if a different velocity field is 
used. If the gas is assumed to be in free-fall toward the 0.5 M Q 
protostar with only thermal line broadening (i.e., no additional 
turbulent velocity field), evidence of a huge abundance jump 
(>1000) can be found for L1448-C. The location of this jump 
is at the 100 K radius of the envelope, where thermal evapora- 
tion can take place. For L1448-C, this radius lies at 33 AU or 
0715, and the interferometer data can be used to test these two 
different interpretations. 

Fig. 9 shows the observed H2CO visibility amplitudes to- 
ward L1448-C. The emission has been averaged over the 
full extent of the line (4-6 kms -1 ). Although the signal-to- 
noise is low, the emission is clearly resolved meaning that hot 
H2CO extends to scales larger than 1". Fig. 9 also presents the 
model predictions assuming a constant para-H 2 CO abundance 
of 6 x 10~ 10 , consistent with the single-dish data (solid line). 
The quality of the fit is measured using a x 2 statistic for those 
visibility amplitudes that are above the zero-expectation limit. 
From the fit to the observed visibilities in Fig. 9 it is evident that 
such a constant H2CO abundance throughout the envelope can- 
not reproduce the interferometer data (x 2 eA = 26.8): the 3o3-2 2 
emission is overproduced on short baselines and underpro- 
duced on the longest baselines, while the 322 — » 22i/3o3 — > 2o2 
model ratio is much lower than the observations. 

One possible explanation is that the para-F^CO abundance 
drastically increases from 6x 10~ 10 in the outer cool parts of the 
envelope to 7 x 10~ 7 when T > 100 K, in accordance with the 
analysis of Maret et al. (2004) and the results for IRAS 16293- 
2422 (Ceccarelli et al. 2000b; Schoier et al. 2002). This situa- 
tion is similar to that for the point source needed to explain the 
continuum emission in §3, since the region where T > 100 K is 
only about 0715 (33 AU) in radius. Enhancing the abundance 
in this region will therefore correspond to adding an unresolved 
point source. The visibility amplitudes for this model are pre- 
sented in Fig. 9 (dashed line) and agree better with observations 
on longer baselines than the constant abundance model, but the 
overall fit is actually slightly worse (x 2 d = 31.4). Allowing for 
a jump at temperatures as low as 50 K cannot be ruled out in 
the case of IRAS 16293-2422 (Schoier et al. 2002; Doty et al. 
2004). This would extend the region of warm material to scales 
of 1" (~ 200) AU, but this still cannot explain the emission on 
the shortest baselines, i.e., on scales of ~ 10 - 20". 

Similarly, the observed emission sampled by the longer 
baselines could be caused, in part, by H2CO emission com- 
ing from the warm circumstellar disk that was suggested to 
explain the compact dust emission in §3.1. Here we adopt a 
disk temperature of 100 K and a corresponding mass (Eq. 1) 
of 0.016 M . Further, assuming a diameter of 70 AU, similar 
to the size of the region where the temperature is higher than 
100 K and consistent with the upper limit from observations, 
the spherically averaged number density of H2 molecules is 
5.3 xlO 9 cirT 3 . While the temperature is similar to that found in 
the inner envelope, the density scale is ~ 10-50 times higher. 
It is found that an H2CO abundance of 4 x 10~ 9 in the disk can 
explain the observed visibilities on the longer baselines. The 
abundance in the disk is about 7 times larger than found from 
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single-dish modelling of the envelope alone, but more than an 
order of magnitude lower than what the thermal evaporation 
model predicts. The H2CO abundance in the disk depends crit- 
ically on the assumed disk properties; e.g., the disk mass used 
above is only a lower limit and the actual value may be an or- 
der of magnitude higher. This would remove the need for a 
H2CO abundance jump all together. Still, since the disk is un- 
resolved this does not alleviate the problem of the relatively 
strong H2CO 322 — » 221 line emission on the shortest baselines. 

In their study of a larger sample of low-mass protostars, 
J0rgensen et al. (2002) show that CO is significantly depleted 
in deeply embedded objects such as L1448-C. They also found 
that intensities of low excitation J = 1 — > lines are not con- 
sistent with constant abundances derived on basis of the higher 
J lines. Similar trends are seen for other molecular species such 
as HCO + and HCN by J0rgensen et al. (2004b), who suggest 
that this is caused by the fact that the time scale for freeze-out 
of CO and other species is longer than the protostellar lifetime 
in the outer envelope. This leads to a 'drop' abundance profile 
(see Fig. 10), with CO frozen out in the cold region of the en- 
velope, but with standard or enhanced abundances in the out- 
ermost low density cloud and in the inner warmer envelope. 
The region over which CO is frozen out is determined by the 
outer radius R2 at which the density is high enough that the 
freeze-out time scale is short compared with the protostellar 
lifetime and the inner radius Ri at which the temperature is low 
enough that CO does not immediately evaporate from the grain 
ice mantles. Photodesorption may also play a role in the outer- 
most region. In order for the freeze-out time scale to be shorter 
than ~ 10 4 years, the density should be higher than ~ 10 5 cirT 3 . 

The H2CO abundance profile is expected to follow at least 
partially that of CO, because destruction of gas-phase CO by 
He + can be a significant source of atomic carbon and oxygen: 

He + + CO -> C + + O + He (2) 

For the typical densities and temperatures in the outer region, 
H2CO is mainly formed through the reaction: 

CH 3 + O -> H 2 CO + H (3) 

so the H2CO abundance should drop in regions where CO is 
frozen out. Indeed, Maret et al. (2004) found a clear correlation 
between the H2CO and CO abundances in the outer envelopes 
where both molecules are depleted, for a sample of eight class 
protostars. Such an effect would show up in the comparison of 
interferometer and single-dish data: the interferometer data are 
mainly sensitive to the 1 - 10" region of the envelope, where 
CO is frozen out. The single-dish lines, however, either probe 
the outer regions of the envelope (low excitation lines), where 
they quickly become optically thick, or the inner regions (high 
excitation lines) which are unaffected by the depletion. 

In order to test this effect several models in which H 2 CO 
is depleted over roughly the same region as CO were consid- 
ered (see Fig. 10). Each model was required to simultaneously 
reproduce all the available multi-transition single-dish data to 
better than the 3cr level. First, an 'anti-jump' model was con- 
sidered where the abundance drops from an initial undepleted 
value, Xq, to X D when the H2 density is larger than 10 5 crrT 3 . 



As can be seen in Fig. 9 using Xq — 5 x 10~ 9 with a drop to 
X D = 6 x 10~ 10 drastically improves the fit to the observed 
3o3 — » 2o2 line emission due to opacity effects. The optically 
thin 322 — » 221 line emission is unaffected by this. The overall 
fit for this model is good, x 2 ed = 2.0. 

Next, a 'drop' H2CO abundance profile was introduced to 
simulate the effects of thermal evaporation in the inner warm 
part of the envelope. First, r ev was taken to be 30 K (see 
discussion in J0rgensen et al. (2002)), roughly the evapora- 
tion temperature of CO. For r ev = 30 K the jump is located 
at Ri = 7.5 x 10 15 cm (Z'3). A para-H 2 CO abundance of 
X Q = 5 x 10~ 10 with a drop to I D = 3 x 10~ 10 in the re- 
gion of CO depletion provides the lowest x 2 and is consistent 
with the single-dish data. However, the fit to the interferome- 
ter data is not good, xLa = 14.6. In particular the 3o3 — > 2o2 
line emission comes out too strong in the model since Xq is 
not allowed to increase enough (constrained by the single-dish 
data) to become optically thick as in the anti-jump model. The 
'drop' model does, however, provide a good description of 
the interferometer emission at both long and short baselines 
if the H2CO abundance remains low out to T > 50 K. For 
T ev = 50 K (R\ = 1.9 x 10 15 cm; 0758), X Q = 5 x 10~ 9 and 
X D = 4 x 10~ 10 provides a good fit with^ ed = 1.4. Raising 
T ev to 100 K (Ri = 5.0 x 10 14 cm; 0715) provides a slightly 
worse fit, xi d = 1.8, for Xq = 4 x 10~ 9 and X D = 4 x 10~ 10 . 
Xq is forced to be in the range 4 - 5 x 10~ 9 from the single 
dish data, so for r ev = 100 K less flux is obtained at shorter 
baselines than compared with the model where r ev = 50 K 
because of the smaller emitting volume. Note that these mod- 
els predict a high 322 — * 22i/3o3 — > 2o2 line ratio at large scales 
where the gas temperature is only ~20 K, due to the fact that the 
3o3 — » 2o2 line in the outer undepleted region becomes optically 
thick. On the longest baselines, compact emission from either 
a disk or an abundance jump would still be consistent with the 
observations. For example, an additional jump of a factor 10 
when T > 100 K (Rq - R\ = 5 x 10 14 cm; X D =4x 10~ 10 ; 
Xq = 4 x 10~ 9 ; Xj = 4 x 10~ 8 ; see Fig. 10) results in an 
additional « 0.1 Jy on all baselines and improves the overal 
fit, xLa - 1 -4> f° r evaporation at this higher temperature (see 
Fig. 9). However, since the observed signal is close to the zero- 
expectation level, no strong conclusions on the presence of this 
additional jump can be made. 

To summarize: while a constant abundance model can ex- 
plain the H2CO emission traced by the single-dish data, it 
underproduces the emission in the interferometer data on the 
longest baselines and produces too much 3o3 -2o2 line emission 
on the shorter baselines. Adding a compact source of emission 
either through a hot region of ice mantle evaporation or a cir- 
cumstellar disk does not provide a better overall fit. A 'drop' 
profile in which the H2CO abundance largely follows that of 
CO alleviates these problems, explaining at the same time the 
emission seen by single-dish and the structure of the emis- 
sion as traced by the interferometer. The relatively high H2CO 
322 - 221 /3o3 - 2o2 ratio at scales of 1 to 10" is in this scenario 
caused by a combination of high optical depth of the 3o3 - 2o2 
line in the outermost region, and a low H2CO abundance in the 
cold dense part of the envelope where CO is frozen out. The 
best-fit abundances in each of these scenarios are summarized 



10 



F. L. Schoier et al.: On the origin of H2CO abundance enhancements in low-mass protostars 




H CO (3 -2 ) 

2 v 03 02 ; 




H CO (3 -2 ) 

2 v 22 21 7 



n 




10 5 -5 

RA offset ["] 



-10 



Fig. 11. OVRO interferometer maps of IRAS 16293-2422 in C ls O 
and H 2 CO line emission (contours) oveiiayed on the 1.37 mm con- 
tinuum emission (greyscale). The molecular line emission has been 
separated into red (dashed lines) and blue (solid lines) components 
(see text for details). Contours start at 1.8 JybeanT 1 kms~' for C ls O 
and at 0.9 Jy beam 1 Ions' 1 for H2CO. The first contour corresponds 
to the '2<x'-level as estimated from the maps and each successive con- 
tour denotes an increase of 2<x. Also indicated is the direction of the 
large scale CO outflow. 



in §5 and compared with those obtained from the similar anal- 
ysis performed for IRAS 16293-2422 in §4.2.2. 

Finally, it should be noted that there may still be other ex- 
planations for the high 322 - 221 /3o3 - 2o2 line ratio on short 
baselines. If the H2CO emission originates from low-velocity 
entrained material in regions where the outflow interacts with 
the envelope, the gas temperature may be increased due to weak 
shocks. Alternatively, the gas temperature can be higher than 
that of the dust due to heating by ultraviolet or X-ray pho- 
tons from the protostar which can escape through the biconical 
outflow cavity and scatter back into the envelope at larger dis- 
tances (cf. Spaans et al. 1995). Detailed quantitative modeling 
of these scenarios requires a good physical model of the heat- 
ing mechanisms and a 2-D radiative transfer and model analy- 
sis, both of which are beyond the scope of this paper. In either 
scenario, the general envelope emission described above still 
has to be added, and may affect the line ratios through opacity 
effects. 

4.2. IRAS 16293-2422 
4.2.1. C 18 emission 

The overall C ls O J = 2 ^ I line emission for IRAS 16293- 
2422 obtained at OVRO is presented in Fig. 11, with the chan- 
nel maps shown in Fig. 12. The emission is clearly resolved 
and shows a ~6" separation between the red (4-7 kms _I ) 
and blue (1-4 kms"') emission peaks. The direction of the 
red-blue asymmetry is roughly perpendicular to the large scale 
CO outflow associated with MM1 (Walker et al. 1988, Mizuno 
et al. 1990, Stark et al. 2004), and may be indicative of overall 
rotation of the circumbinary material encompassing both MM1 
and MM2. The morphology and velocity structure is also con- 
sistent with a large (~1000 AU diameter) rotating gaseous disk 
around just MM1, however. Such large rotating gaseous disks 
have been inferred around other sources, including the class 
I object L1489 (Hogerheijde 2001) and the classical T Tauri 
star DM Tau (Dutrey et al. 1997), although their C ls O emis- 
sion is usually too weak to be detected. The C 18 emission 
clearly avoids the region between the binaries, consistent with 
the conclusion from the continuum data that this region is void 
(see §3.2). 

In Fig. 13 the interferometer data, restored with a 22" 
beam, are compared with the JCMT single-dish flux. The in- 
terferometer only recovers ~5% of the total single-dish flux, 
mainly at extreme velocities. The interferometer is not sensi- 
tive to the large scale static emission close to the cloud velocity 
of «4 kms" 1 . 

The C 18 data can be used to further test the density struc- 
ture of the circumbinary envelope. However, as discussed in 
§4. 1 there is now growing evidence that CO is depleted in the 
cool outer parts of the envelope so that extra care is needed 
when deducing the density structure from CO observations 
alone. Schoier et al. (2002) were able to reproduce single- 
dish CO observations using a constant abundance of about 
3 + 1 x 10~ 5 throughout the envelope and with jump mod- 
els at 20 K, i.e., the characteristic evaporation temperature, 
Tqo, of pure CO ice. Recently, Doty et al. (2004) found that 
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Fig. 12. OVRO interferometer maps of C' s O J = 2 — > 1 line emission (contours) from IRAS 16293-2422. The C 18 emission has been 
separated into bins of 1 kms 1 and contours are in steps of 0.8 JybeanT 1 kms~' (2a). Dotted contours indicate negative values. In the panel 
with the velocity-integrated line intensity the contours start at 2.0 Jy beam -1 km s -1 (2cr). Also indicated are the positions of the two continuum 
sources as well as the red- (dashed) and blue-shifted (solid) parts of the large scale outflow. 



7co ~ 20 K provided the best fit in their chemical modelling of 
a large number of molecular species. Here the C ls O J — 2 — > 1 
data are analyzed assuming that the emission originates from 
1) a 'standard' envelope centered on MM1; and 2) a circumbi- 
nary envelope with a cavity centered between the positions of 
the protostars MM1 and MM2. 

Fig. 14 shows the result of applying the (u, v) sampling 
of the observations to the C O envelope models. In both the 
model centered on MM1 and the common envelope model 
(with a cavity) the observed visibilities are relatively well re- 
produced using the same constant C 18 abundance of 6 x 1CT 8 
(solid lines) as obtained from the single-dish analysis, except 
perhaps at the longest baselines. Note that the density and tem- 
perature structures are slightly different for the model with 
a cavity (see §3.2). Next, thermal evaporation models with a 
drastic jump at Tco = 20 K, as suggested by Doty et al. (2004), 
were considered. The abundance in the region of depletion, 
i.e., when T < 20 K and «h 2 > 1 x 10 5 crrT 3 , was fixed to 
2 x 10~ 8 , the upper limit obtained by Doty et al. (2004). An 



undepleted C 18 abundance, Xq, of 8 x 10~ 8 (corresponding to 
a total CO abundance of about 5 x 10~ 5 ) is found to provide 
reasonable fits to the observed visibilities in Fig. 14 (dashed 
lines) as well as CO single-dish data, and is also just consistent 
with the chemical modelling performed by Doty et al. (2004). 
This CO abundance is a factor of 2 to 5 lower than typical un- 
depleted abundances. Assuming a still higher evaporation tem- 
perature of 50 K brings X for C ls O up to 2 x 10~ 7 (CO up to 
1 . 1 x 10 -4 ), more in line with what is typically observed in dark 
clouds. Doty et al. (2004) could not rule out such high evapo- 
ration temperatures in their chemical modelling. In all, the suc- 
cess of the envelope model in explaining both the continuum 
emission as well as the observed CO emission is encouraging 
and will aid the interpretation of the H2CO observations. 

The flux detected on the longer baselines provides a limit 
on the maximum CO abundance in the disk. For a disk around 
MM1 with diameter 250 AU, temperature 100 K and the mass 
0.09 M Q (see §3.2) the spherically averaged H2 number den- 
sity is 6.4 x 10 8 cm 4 . It is found that a C I8 abundance of 
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Fig. 13. Comparison between the C I8 J = 2 — > 1 line emission to- 
wards IRAS 16293-2422 at the source position, from JCMT single- 
dish observations (line diagram) and OVRO interferometric observa- 
tions (histogram) restored with the JCMT beam (22"). The OVRO 
spectrum has been scaled in order to account for the flux seen in the 
JCMT spectrum. 



w 2 x 10~ 7 , corresponding to a total CO abundance of about 
1.1 x 10~ 4 , can account for the emission on the longest base- 
lines. Assuming a lower temperature in the disk of 40 K and 
the correspondingly higher mass of 0.24 M gives almost the 
same upper limit to the amount of CO in the disk. 

4.2.2. H 2 CO emission 

The single-dish observations (see Fig. 15) show that the H2CO 
line emission is extended to scales of ~30". The single-dish 
322 — > 22i/3o3 — > 2o2 line ratio is »0.2, suggesting that a cold 
(30-40 K) envelope component dominates the single-dish flux. 
In contrast, the interferometer 322 — * 22i/3(>3 — » 2o2 line ra- 
tio is 0.69 ± 0.23 for the red-shifted emission near MM1 and 
0.75 + 0.32 for the blue-shifted emission close to MM2. This 
indicates that the temperature is in excess of ~ 150 K assum- 
ing the density to be at least 10 6 cm 3 . However, as noted for 
L1448-C, optical depth effects and abundance gradients can 
affect this ratio. 

The velocity channel maps of the H2CO 3o3 — * 2o2 interfer- 
ometer line emission obtained at OVRO are shown in Fig. 16 
for IRAS 16293-2422, whereas the total H 2 CO 3 22 -> 2 2 i 
and 3rj3 — > 2q2 maps are included in Fig. 11. The emission 
is clearly resolved and shows a ~6" separation between the red 
(4-7 kms -1 ) and blue (1-4 kms -1 ) emission peaks. The di- 
rection of the red-blue asymmetry is again roughly perpendic- 
ular to the large scale CO outflow associated with MM1, but in 
contrast with the C 18 OVRO maps, no blue-shifted emission 
close to MM1 is observed. Instead, the strongest blue-shifted 
emission occurs close (<1") to MM2. The red-shifted H2CO 
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Fig. 14. Visibility amplitudes of the observed C ls O 2 — > 1 line emis- 
sion obtained at OVRO toward IRAS 16293-2422 as functions of 
the projected baseline length from the phase center, taken to be at 
(2", -3") for the envelope around MM1 (top) and (072, -1'.' 2) for 
the common envelope (bottom). The observations, averaged over 1 
to 7 km s and binned to 10 kA, are plotted as filled symbols with lcr 
error bars. The dotted histogram represents the zero-expectation level. 
Also shown is the result of applying the same (u, v) sampling to the cir- 
cumstellar model using a constant C 18 abundance of 6 x 10~ 8 (solid 
lines) for an envelope centered on the protostar MM1 and a common 
envelope with a cavity. Envelope models where CO is depleted in a 
region where nn 2 > 1 x 10 5 cirT 3 and T < T ev are also indicated as 
dashed (T ev = 20 K; X = 8 x 10~ 8 ; X D = 2 x 10~ 8 ) and dot-dashed 
(r cv = 50 K; X = 2 x 10~ 7 ; X D = 2 x 10~ 8 ) lines. 

emission is found to the south of MM1, but again peaks closer 
to MM1 than the C I8 emission does. 

As for C ls O, only a small fraction (~ 5-20%) of the single- 
dish flux is recovered by the interferometer (see Fig. 17). The 
shapes of the H2CO 3(>3 — > 2o2 and 322 — > 221 lines are very 
similar to that of C ls O J = 2 — > 1 (see Fig. 13), when the 
interferometer data are restored with the JCMT beam. 

The H2CO emission is interpreted in terms of a common 
envelope scenario using the cavity model that was seen to best 
reproduce the emission in §3.2. Alternative models with the 
envelope centered on MM1 or MM2 have been considered as 
well, but lead to the same overall conclusions. H2CO models 
with and without any abundance jumps or drops are produced 
for this physical model, following the analysis performed for 
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Fig. 15. JCMT single-dish spectral maps of the H 2 CO emission toward 
IRAS 16293-2422. The velocity scale is in kms" 1 and the intensity is 
the main beam brightness temperature in K as indicated in the upper 
right subplot. The spectral resolution is 0.22 kms -1 . The 3 22 — * 2 2 i 
line data have been smoothed to a two times lower spectral resolution. 



L1448-C (see Fig. 10). Applying the (u, v) sampling from the 
observations to these models produces visibility amplitudes 
that are compared with observations in Fig. 18. 

Using a constant para-H 2 CO abundance of 5 x 10~ 10 de- 
rived from the single-dish modelling performed in Schoier et al. 
(2002) produces too much flux on the shorter baselines for both 
the 3 03 -> 2 02 and 3 22 -» 2 2 i transitions. For IRAS 16293- 
2422, an abundance jump in the inner hot core due to evapo- 
ration of ice mantles is well established from multi-transition 
single-dish modelling (van Dishoeck et al. 1995; Ceccarelli 



et al. 2000b; Schoier et al. 2002). Schoier et al. (2002) con- 
strain the location of the jump to > 40 K, with a jump in abun- 
dance of one to three orders of magnitude depending on the 
adopted evaporation temperature. It is found that a jump model 
with T ev = 50 K (at /?! = 6.2 x 10 15 cm = 410 AU = 276) 
improves the fit to the data (xL* = 7.1) using Xq = 1 x 10~ 10 
and X s = 3 x 10~ 9 (Fig. 18, left panels). However, this model 
still produces too much 3o3 — > 2 ( )2 line emission on the short- 
est baselines. This is similar to L1448-C where a better fit was 
obtained by letting the abundance increase again in the outer 
region when n H2 < 1 x 10 5 cm -3 to provide significant opti- 
cal depth in the 3o3 — > 2o2 transition. Such an anti-jump model 
with X D = 1 x 10~ 10 and X Q = 3 x 10~ 9 also improves the fit 
(xl d = 3.6) in the case of IRAS 16293-2422. 

In the drop models two different evaporation temperatures 
are considered: r ev = 50 K and 30 K (at R\ = 1.1 x 10 16 cm 
= 740 AU = 476). As shown in Fig. 18 (right panels) the 30 K 
model does not provide a good fit to the data (xLa - 17.5). 
Using r ev = 50 K instead allows Xq to be higher and a near- 
perfect fit (x 2 mi = 0.4) can be found to both the 3o3 -> 2q 2 
and 322 —> 2n line emission. In this case Xu = 1 x 10~ 10 and 
Xo = 4x 10~ 9 , i.e., the jump is of a factor 40. A similar jump at 
50 K was found by Ceccarelli et al. (2001) from their analysis 
of the H2CO single-dish data. The best-fit abundances in each 
of these scenarios are summarized in §5 and compared with 
those obtained from the similar analysis performed for L1448- 
C in §4.1. 

The observed flux at the longest baselines gives an up- 
per limit to the para-F^CO abundance in a disk around MM1. 
Using the properties of the disk as in §4.2.1 an abundance of 
1 x 10 -9 is found to produce about 0.5 Jy on all baselines. 
For the «30% more massive disk around MM2 a slightly lower 
value for the para-F^CO abundance is found. 

Given the complexity of the H2CO emitting region on small 
scales (<10") the envelope model presented in Schoier et al. 
(2002) is not adequate to fully describe the morphology of the 
emission observed by the interferometer. Here we simply note 
that the envelope model can explain the observed flux if jumps 
in abundance are introduced. The fact that both the H2CO and 
C 18 observations can be well reproduced by the common en- 
velope model with a cavity further supports the finding from 
the continuum observations that little material seems to be lo- 
cated at inter-binary scales. 

5. Origin of the H 2 CO emission 

Here the results from the previous sections are summarized. 
Various competing scenarios on the origin of the observed 
H2CO emission are discussed and predictions for future gen- 
eration telescopes are presented, illustrating their potential to 
distinguish competing scenarios. 

5.1. Envelope and/or outflow emission? 

In §4 the physical envelope models derived from single-dish 
observations, and tested against interferometric continuum ob- 
servations in §3, have used as a basis for interpreting the ob- 
served H 2 CO emission. The results are summarized in Table 2. 
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Fig. 16. OVRO interferometer maps of H2CO 3o3 — » 2o2 line emission (contours) for IRAS 16293-2422. The H2CO emission has been separated 
into bins of 1 kms 1 and contours are in steps of 0.4 JybeairT 1 kms" 1 {2a). In the panel with the velocity-integrated line intensity the contours 
start at 1.0 Jy beam" 1 kms" 1 (2<x). Also indicated are the positions of the two continuum sources as well as the red- (dashed) and blue-shifted 
(solid) parts of the outflow. 



It is found that for both IRAS 16293-2422 and L1448-C, the 
best fit to the H2CO interferometric observations is obtained 
with a 'drop' abundance profile, in which the H2CO abundance 
is lower by more than an order of magnitude in the cold dense 
zone of the envelope but is high in the inner- and outermost re- 
gions. The outer radius of this 'drop' -zone is set by the distance 
at which the density drops below 10 5 cm" 3 ; the inner radius by 
the distance at which the temperature is above the evaporation 
temperature r ev . Indeed, such a 'drop' model with r ev =50 K 
gives very good x 1 fits and reproduces both the single-dish 
and interferometer data. The fact that two lines of the same 
molecule with different optical depths were observed simul- 
taneously with OVRO was essential to reach this conclusion. 
The presence of abundance enhancements in regions where 
T > 50 K is consistent with the detailed analysis of multi- 
transition single-dish H2CO observations, although the actual 
values derived here are somewhat different. Interestingly, the 
inferred abundances X D and Xq for the best-fit 50 K drop mod- 
els are comparable for the two sources. The presence of addi- 



tional jumps with Xj > Xq in the innermost hot core region or 
disk cannot be established with the current data, but requires in- 
terferometer observations of higher excitation lines (see §5.4). 

Can some of the enhanced H2CO originate in the outflow? 
The red-blue asymmetry observed for L1448-C is consistent 
with the high velocity outflow seen in CO and SiO, but the 
velocities seen for H2CO are significantly lower. One explana- 
tion is that the emission originates in the acceleration region 
of the outflow, estimated to be within 2" radius from the star 
(Guilloteau et al. 1992). An alternative, more plausible scenario 
is that if H2CO is associated with the outflow, but located in 
low-velocity entrained material in regions where the outflow 
interacts with the envelope, since the red-shifted emission ap- 
pears to be extended to scales much larger than 2". This would 
be similar to the case of HCO + (Guilloteau et al. 1992). 

For IRAS 16293-2422 the morphology indicates emission 
in a rotating envelope perpendicular to the direction of the 
(large-scale) outflow. However, the complicated physical struc- 
ture of this protobinary object on small spatial scales is not 
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Fig. 18. Visibility amplitudes of the observed H2CO line emission obtained at OVRO toward IRAS 16293-2422 as functions of the projected 
baseline length from the phase center, taken to be at (072, -1'.'2) for the common envelope. The observations, averaged over 1 to 7 km s" 1 and 
binned to 10 kA, are plotted as filled symbols with lcr error bars. The dotted histogram represents the zero-expectation level. Also shown are 
the results of applying the same (u,v) sampling to the circumbinary envelope model using various H2CO abundance distributions. All models 
are consistent with available multi-transition single-dish data to within 3cr. See text and Fig. 10 for details on these models. 



Table 2. Best fit H 2 CO envelope models". 
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" All abundances refer to para-H 2 CO only. 

* Reduced^ 2 from interferometer data. All models are consistent with available single-dish data to within the 3<x level. 
c Jump at 100 K for L1448-C and 50 K for IRAS 16293-2422. 



well represented by our spherically symmetric model, so it is 
not possible to rule out a scenario where some of the emission 
originates on larger scales due to interaction with the outflow. 

To quantify the role of outflows in producing H2CO abun- 
dance enhancements and liberating ice mantles, H2CO interfer- 
ometer data at ~1" resolution for a larger sample of sources are 
needed to investigate whether the velocity pattern is system- 
atically oriented along the outflow axis as in L1448-C. Also, 
higher sensitivity could reveal whether the profiles have more 
extended line wings. 



5.2. Photon heating of the envelope? 

The current models assume that the gas temperature equals the 
dust temperature. Detailed models of the heating and cooling 
balance of the gas have indicated that this is generally a good 
assumption within a spherically symmetric model (Ceccarelli 
et al. 1996; Doty & Neufeld 1997). However, if the gas tem- 
perature were higher than the dust temperature in certain re- 
gions, this would be an alternative explanation for the increased 
322 - 22i/3o3 - 2o2 ratio in the interferometer data on short base- 
lines. One possibility discussed in §4.1 would be gas heating 
by ultraviolet (UV)- or X-ray photons which impact the in- 



16 



F. L. Schoier et al.: On the origin of ITCO abundance enhancements in low-mass protostars 



2.0 



1.5 



0.5 



0.0 



-0.5 



10 



IRAS 16293-2422 
H CO (3 -2 ) 

2 v 22 21 y 



OVROxlO 




10 



15 



t4 



IRAS 16293-2422 
H CO (3 -2 ) 

2 v 03 02' 



0VR0X15 




5 10 

Velocity [km s -1 ] 



15 



Fig. 17. Comparison between the H2CO line emission towards IRAS 
16293-2422 at the source position, from JCMT single-dish obser- 
vations (line diagram) and OVRO interferometric observations (his- 
togram) restored with the JCMT beam (22"). The OVRO spectrum 
has been scaled in order to account for the flux seen in the JCMT 
spectrum. 



ner envelope and can escape through the biconical cavity ex- 
cavated by the outflow. If such photons scatter back into the 
envelope, they can raise the gas temperature to values signifi- 
cantly in excess of the dust temperature in part of the outer en- 
velope (e.g., Spaans et al. 1995). For IRAS 16293-2422, such 
photons can further escape through the circumbinary cavity, so 
that the photons from e.g., MM2 can influence the inner enve- 
lope rim around MM1. Since this model would affect the ex- 
citation of all molecules present in this gas, not only H2CO, 
it can be tested with future multi-line interferometer data of 
other species. Moreover, the presence of UV photons would 
have chemical consequences producing enhanced abundances 
of species like CN in the UV-affected regions, which should be 
observable. Note that in this scenario, the general colder enve- 
lope still has to be added, which, as noted above, can affect the 
line ratios. 



5.3. Disk emission? 

The detailed modelling of the continuum emission performed 
in §3 reveals that there is compact emission in both IRAS 
16293-2422 and L1448-C that cannot be explained by the en- 
velope model. The most likely interpretation is that of accretion 
disks. For L1448-C, where the inner region appears to be less 
complex, it is found that the observed compact H2CO emis- 
sion can be equally well explained originating from a disk as 
from the inner hot core. However, the need for a drastic jump 
in abundance depends critically on the properties of the disk. 
An upper limit on the para-H2CO abundance in the disk of 
~4xl0~ 9 is derived adopting a temperature of 100 K, a mass of 
0.016 M and a size of 70 AU and assuming that all the flux on 
long baselines arises is due to the disk. For IRAS 16293-2422 
an upper limit of the abundance in the MM1 disk of ~lxl0~ 9 
is obtained using a mass of 0.09 M Q and a size of 250 AU. A 
slightly lower upper limit is obtained for MM2 using a disk 
mass of 0.12 M . Moreover, the difference between the mor- 
phology of the C ls O and H2CO emission seems to indicate 
that at least the disk around MM1 contains little H2CO, even 
though CO may be nearly undepleted. Note that geometrical 
effects, for example from the disk potentially shielding parts of 
the envelope, can cause differences in the appearance between 
MM1 and MM2. 

H2CO has been detected in protoplanetary disks of T Tauri 
stars (Dutrey et al. 1997; Aikawa et al. 2003), where confu- 
sion due to an envelope or outflow is negligible. Thi et al. 
(2004) find an integrated H2CO 3o3 — > 2(>2 line intensity of 
0.14 Kkms 1 (about 1.3 Jykms 1 ) for the T Tauri (class II) 
star LkCa 15 using the IRAM 30 m telescope. The beam- 
averaged (1078) H 2 CO abundance is about 4 x 10" 11 . The 
OVRO data for IRAS 16293-2422 presented here give a fac- 
tor of about 50 stronger emission. For L1448-C the emission 
is about a factor of 10 stronger after correcting for the distance. 
Thus, if the compact H2CO emission were coming from disks, 
they would have to be 'hotter' or have higher abundances than 
in the class II phase. Accretion shocks in the early stages could 
be responsible for such increased disk temperatures. High- 
angular resolution (< 1") data on the velocity structure of the 
H2CO lines are needed to distinguish the disk emission from 
that of the inner envelope. 

5.4. Predictions for future generation telescopes 

In Table 3, the predicted H2CO line intensities, integrated over 
the full extent of the line, are presented for L1448-C using: 1) 
the envelope model with a constant abundance of 6 x 10~ 10 , 
2) introducing a jump when T > 100 K (Xj = 1.5 X 10~ 7 ; 
X D = 6 x 10~ 10 ), 3) a 'drop' abundance profile where H 2 CO 
is depleted over the same region as CO (r ev = 50 K; Xq = 
5 x 10~ 9 ; X D = 4 x 10~ 1() ); and 4) the envelope + disk model 
(X D = 6 x 10~ 10 and disk parameters from §4.1). Beamsizes 
of 073 and 3" were assumed to characterize the typical spa- 
tial resolutions of current and future interferometers. The cor- 
responding intensities picked up by a single-dish beam of 15" 
are shown for comparison. Using single-dish data alone, it will 
be difficult to discriminate between these competing scenarios 
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" Energy of the upper energy level involved in the transition. 

* Constant para-H 2 CO abundance of 6 x 10~ 10 throughout the envelope. 

A jump in abundance of a factor 250, from 6 x 10~ 10 , when T > 100 K . 
d H2CO is depleted by about an order of magnitude over roughly the same region as CO (r c v = 50 K; Xq = 5 x 10~ 9 ; Xv = 4.0 x 10~ 10 ) 
e H 2 CO abundance in the disk is 25 times higher than that in the envelope of 6 x 10~ 10 . 



unless the highest frequency lines are obtained. Current inter- 
ferometers such as OVRO working in the 1 mm window can 
however constrain some characteristics of the abundance varia- 
tions in the envelopes, such as the presence of a drop abundance 
profile. This seems to be the case for the two sources studied 
here, IRAS 16293-2422 and L1448-C. 

Finally, it is clear that observations at 0."3 will have the po- 
tential to further discriminate a hot core scenario from that of 
a warm disk. Indeed, the much improved resolution and sen- 
sitivity of next generation interferometers such as CARMA, 
(e)-SMA, (upgraded) IRAM and ALMA will greatly aid in dis- 
tinguishing between the competing scenarios discussed above. 
They will also provide additional constraints on the morphol- 
ogy and velocity structure of the gas on larger and smaller 
scales. In addition to searching for outflow motions, it is of 
considerable interest to determine if the hot core gas shows any 
evidence of infall motions toward the cental source(s). If chem- 
ically processed material such as H2CO and other organics is in 
a state of infall toward the central protostars it will likely be in- 
corporated in the growing protoplanetary disk and become part 
of the material from which planetary bodies are formed. With 
the present data, it is not possible too uniquely separate infall 
motions from those of rotation and/or outflow. 

6. Conclusions 

A detailed analysis of the small scale structure of the two 
low-mass protostars IRAS 16293-2422 and L1448-C has been 
carried out. Interferometric continuum observations indicate 
that the inner part of the circumbinary envelope around IRAS 
16293-2422 is relatively void of material on scales smaller 
than the binary separation (~5")- This implies that the clear- 
ing occurs at an early stage of binary evolution and that IRAS 
16293-2422 may well develop into a GG Tau-like object in the 
future. The bulk of the observed emission for both sources is 
well described using model envelope parameters constrained 
from single-dish observations, together with unresolved point 
source emission, presumably due to circumstellar disk(s). 

Simultaneous H2CO line observations indicate the presence 
of hot and dense gas close to the peak positions of the contin- 
uum emission. For both IRAS 16293-2422 and L1448-C, the 



observed emission cannot be reproduced with a constant abun- 
dance throughout the envelope. The H2CO 322 - 221 /3o3 - 2o2 
ratio on short baselines (2 - 10") is best fit for both sources by 
an H2CO 'drop' abundance profile in which H2CO, like CO, 
is depleted (by more than an order of magnitude) in the cold 
dense region of the envelope where T < 50 K, but is relatively 
undepleted in the outermost region where « H2 < 1 x 10 5 cirT 3 . 
In the inner region for T > 50 K, the abundance jumps back 
to a high value comparable to that found in the outermost un- 
depleted part. Additional H2CO abundance jumps — either in 
the innermost 'hot core' region or in the compact circumstellar 
disk — cannot be firmly established from the current data set. 

Based on the morphology and line widths, little of the ob- 
served emission toward IRAS 16293-2422 is thought to be di- 
rectly associated with the known outflow(s). Instead, the emis- 
sion seems to be tracing gas in a rotating disk perpendicular to 
the large scale outflow. For L1448-C, the morphology of the 
H2CO line emission is consistent with the outflow, but the line 
widths are significantly smaller and the emission is extended 
over a large area. Although the above envelope model with a 
'drop' abundance profile can fit the observations, other scenar- 
ios cannot be ruled out. These include the possibility that the 
outflow (weakly) interacts with the envelope producing regions 
of enhanced density and temperature in which some H2CO is 
liberated and entrained, and a scenario in which part of the en- 
velope gas is heated by UV- or X-ray photons escaping through 
the outflow cavities. 

It is clear that high sensitivity and high spatial resolution 
observations are critically needed for a better understanding 
of the chemical and dynamical processes operating during star 
formation. At the same time, this study also demonstrates that 
detailed radiative transfer tools are essential for a correct anal- 
ysis of such data. In the future, the modeling codes need to be 
extended to multiple dimensions (>1D) in order to fully tackle 
the complex geometry hinted at in current data sets. Predictions 
for future arrays are provided that illustrate their potential to 
discriminate between competing scenarios for the origin of the 
H2CO abundance enhancements, and presumably that of other 
complex organics, in low-mass protostars. 
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